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We report on a search for inclusive production of squarks and gluinos in pp collisions at ^/s — 
1.96 TeV, in events with large missing transverse energy and multiple jets of hadrons in the final 
state. The study uses a CDF Run II data sample corresponding to 2 fb _1 of integrated luminosity. 
The data are in good agreement with the standard model predictions, giving no evidence for any 
squark or gluino component. In an R-parity conserving minimal supergravity scenario with Ao = 0, 
pL < and tan/3 = 5, 95% C.L. upper limits on the production cross sections in the range between 
0.1 pb and 1 pb are obtained, depending on the squark and gluino masses considered. For gluino 
masses below 280 GeV/c 2 , arbitrarily large squark masses are excluded at the 95% C.L., while for 
mass degenerate gluinos and squarks, masses below 392 GeV/c 2 are excluded at the 95% C.L. 
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Supersymmetry (SUSY) Q 

is regarded as a possible 
extension of the standard model (SM) that naturally 
solves the hierarchy problem and provides a possible can- 
didate for dark matter in the Universe. SUSY introduces 
a new symmetry that relates fermionic and bosonic de- 
grees of freedom, and doubles the SM spectrum of parti- 
cles by introducing a new supersymmetric partner (spar- 
ticle) for each particle in the SM. Results on similar in- 
clusive searches for SUSY using Tevatron data have been 
previously reported by both the CDF and DO experi- 
ments in Run I 2] and by the DO experiment in Run 
II Q. This Letter presents new results on an inclusive 
search for squarks and gluinos, supersymmetric partners 
of quarks and gluons, based on data collected by the CDF 
experiment in Run II and corresponding to 2.0 fb _1 of 
integrated luminosity. The analysis is performed within 
the framework of minimal supergravity (mSUGRA) Q 
and assumes R-parity conservation where sparticles are 
produced in pairs and the lightest supersymmetric parti- 
cle (LSP) is stable, neutral, and weakly interacting. The 
expected signal is characterized by the production of mul- 
tiple jets of hadrons from the cascade decays of squarks 
and gluinos and large missing transverse energy lfi T [5| 
from the presence of two LSPs in the final state. In a 
scenario with squark masses M~ significantly larger than 
the gluino mass M~ at least four jets in the final state are 
expected, while for M~ > M~ dijet configurations domi- 
nate. Separate analyses are carried out for events with at 
least two, three, and four jets in the final state and with 
different requirements on the minimum ]fi T . The results 
are compared to SM background predictions from quan- 
tum chromodynamics (QCD) jets, W and Z/~f* bosons 
with accompanying jets (VF+jets and Z/7*+jets), top 
quark, and diboson (WW, ZW and ZZ) processes. 

The CDF II detector is described in detail else- 
where @. The detector has a charged particle tracking 
system that is immersed in a 1.4 T solenoidal magnetic 
field coaxial with the beam line, and provides coverage in 
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the pseudorapidity [B[ range |?7| < 2. Segmented sampling 
calorimeters, arranged in a projective tower geometry, 
surround the tracking system and measure the energy 
of interacting particles for < 3.6. The central elec- 
tromagnetic and hadronic calorimeters cover the region 
\rj\ < 1, while the end wall hadronic calorimeter provides 
coverage out to \r]\ < 1.3. Forward electromagnetic and 
hadronic calorimeters cover the regions 1.1 < \rj\ < 3.6 
and 1.3 < \rj\ < 3.6, respectively. Cherenkov counters 
in the region 3.7 < |?/| < 4.7 measure the number of 
inelastic pp collisions to determine the luminosity [7[ • 

Simulated event samples are used to determine de- 
tector acceptance and reconstruction efficiency, estimate 
SM background contributions, and compute the number 
of expected SUSY signal events. Samples of simulated 
QCD-jets, ti production, and diboson processes are gen- 
erated using the PYTHIA 6.216 [8|] Monte Carlo genera- 
tor with Tune A The normalization of the QCD- 
jets sample is extracted from data in a low Ifi T region, 
while ti and diboson samples are normalized to next- 
to-leading order (NLO) predictions 1(J HI]- Samples of 
simulated Z/7*+jets and iy+jets events are generated 
using the ALPGEN 2.1 program [13] where exclusive sub- 
samples with different jet multiplicities are combined, 
and the resulting samples are normalized to the mea- 
sured Z and W inclusive cross sections [13j . Finally, sam- 
ples of single topevents are produced using the MADE- 
VENT program [14| and normalized using NLO predic- 
tions [15j]. In mSUGRA, the mass spectrum of sparticles 
is determined by five parameters: the common scalar 
and gaugino masses at the GUT scale, Mq and Mjy 2 ) 
respectively; the common trilinear coupling at the GUT 
scale, Aq; the sign of the Higgsino mixing parameter, 
fj,; and the ratio of the Higgs vacuum expectation val- 
ues, tan/3. The mSUGRA samples are generated using 
the ISASUGRA implementation in PYTHIA with A = 0, 
fi < 0, and tan/3 = 5, as inspired by previous studies [lij ]. 
A total of 132 different squark and gluino masses are 
generated via variations of Mq and M1/2 in the range 
Ma < 600 GeV/c 2 and 50 < M 1/2 < 220 GeV/c 2 . At 
low tan/3, the squarks from the first two generations are 
nearly degenerate, whereas the mixing of the third gener- 
ation leads to slightly lighter sbottom masses and much 
lighter stop masses. In this analysis, stop pair produc- 
tion processes are not considered. The contribution from 
hard processes involving sbottom production is almost 
negligible, and is not included in the calculation of the 
signal efficiencies to avoid a dependency on the details 
of the model for squark mixing. The mSUGRA samples 
are normalized using NLO cross sections as determined 
by PROSPINO 2.0 [l7j, with input parameters provided by 
isajet 7.74 fill CTEQ61M parton distribution func- 
tions (PDFs) [l9( are used, and renormalization and fac- 
torization scales are set to the average mass of the 
sparticles produced in the hard interaction. The Monte 
Carlo events are passed through a full CDF II detector 
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simulation (based on GEANT3 [2l| and GFLASH [24]) and 
reconstructed and analyzed with the same analysis chain 
as for the data. 

Data are collected using a three-level trigger system 
that selects events with $J T > 35 GeV and at least two 
calorimeter clusters with Et above 10 GeV. The events 
are then required to have a primary vertex with a z po- 
sition within 60 cm of the nominal interaction. Jets are 
reconstructed from the energy deposits in the calorime- 
ter towers using a cone-based jet algorithm [23[ with 
cone radius R = y A77 2 + A</> 2 = 0.7, and the measured 
is corrected for detector effects and contributions 
from multiple pp interactions per crossing at high in- 
stantaneous luminosity, as discussed in Ref. [24|. The 
events are required to have at least two, three, or four 
jets (depending on the final state considered), each jet 
with corrected transverse energy E^ > 25 GeV and 
pseudorapidity in the range |^ ct | < 2.0, and at least 
one of the jets is required to have |?yj ct | < 1.1. Fi- 
nally, the events are required to have ~$ T > 70 GeV. 
For the kinematic range in ~$ T and the of the 

jets considered in this analysis, the trigger selection is 
100% efficient. Beam- related backgrounds and cosmic 
rays are removed by requiring an average jet electro- 
magnetic fraction /„ = Ejets E T*em/ Ejets > 0.15, 
where E^ em denotes the electromagnetic component of 
the jet transverse energy, and the sums run over all the 
selected jets in the event. In addition, the events are 
required to have an average charged particle fraction 
foh = EjetePT.Vk/Ejete-^T* > 0-15, where, for each se- 
lected jet with 1 77^' I < 1.1, Pr\rk ^ s computed as the 
scalar sum of the transverse momenta p^ ack of tracks 
with p£ ack > 0.3 GeV/c and within a cone of radius 
R = 0.4 around the jet axis. The requirements on / em 
and / c h reject events with anomalous energy deposition 
in the hadronic section of the calorimeter or energy de- 
posits in the calorimeter inconsistent with the observed 
activity in the tracking system, and have no significant 
effect on the mSUGRA and SM Monte Carlo samples. 

The dominant QCD-jets background with large $ T 
originates from the misreconstruction of the jet ener- 
gies in the calorimeters. In such events the $ T di- 
rection tends to be aligned, in the transverse plane, 
with one of the leading jets in the event. This back- 
ground contribution is suppressed by requiring an az- 
imuthal separation A</>($ T — jet) > 0.7 for each of the 
selected jets in the event. In the case of the four-jets 
analysis, the requirement for the least energetic jet is 
limited to A0(^ T — jet) > 0.3. Finally, in the two- 
jets analysis case the events are rejected if they con- 
tain a third jet with E^ > 25 GeV, |?7 jot | < 2.0, and 
/S.<p{$ T — jet) < 0.2. The SM background contributions 
with energetic electrons [25| in the final state from Z and 
W decays are suppressed by requiring E^ em /E^ < 0.9 
for each selected jet in the event. In addition, events 



that have one isolated track with p^ ack > 10 GeV/c and 
A<fi(]£ T — track) < 0.7, or two isolated tracks with an 
invariant mass 76 < < 106 GeV/c 2 , are vetoed to 

reject backgrounds with W or Z bosons decaying into 
muon or tau leptons. 

An optimization is carried out to determine, for each 
final state, the lower thresholds on the E^ of the 

individual jets, and Ht, defined as Ht = Ejets^T* 
For each mSUGRA sample, the procedure maximizes 
S/ VB, where 5* denotes the number of SUSY events and 
B is the total SM background. The results from the dif- 
ferent mSUGRA samples are then combined to define, 
for each final state, a single set of lower thresholds that 
maximizes the search sensitivity in the widest range of 
squark and gluino masses (see Table I). As an example, 
for M~ = M~ and masses between 300 GeV/c 2 and 

400 GeV/c 2 , values for S/VB in the range between 20 
and 6 are obtained, corresponding to SUSY selection ef- 
ficiencies of 4% to 12%, respectively. 



Lower thresholds (GeV) 



Final state $ T 



77 rfct(l) p jct(2) pjct(3) p jct(4) 

&T+ > 2 jets 180 330 165 100 

i T + > 3 i cts 120 330 140 100 25 

£ T + > 4 jets 90 280 95 55 55 25 



TABLE I: Optimized lower thresholds on lp T , Ht, and 
Ei^ (i = 1 — 4) for each analysis. 



A number of control samples in data are considered to 
test the validity of the SM background predictions, as ex- 
tracted from simulated events. The samples are defined 
by reversing the logic of some of the selection criteria 
described above. A sample dominated by QCD jets is 
obtained by requiring that at least one of the selected 
jets is aligned with the ]fi T direction. A control sample 
dominated by Z/7*+jets, VU+jets, and top-quark pro- 
cesses with highly energetic electrons in the final state 
is obtained after requiring 

^Tom/^T* > °- 9 for at least 
one of the jets. Similarly, a control sample with highly 

energetic muons in the final state is created by requiring 

the presence of an isolated track with p^, ack > 10 GeV/c 

and Acj)(]fi T — track) < 0.7, or two isolated tracks with 

76 < M tr ks < 106 GeV/c 2 . Good agreement is observed 

between the data and the SM predictions in each of the 

control regions for all the final states considered. 

A detailed study of the systematic uncertainties is car- 



ried out for each final state [27J, |28|. A 3% uncertainty 
on the absolute jet energy scale [24| in the calorime- 
ter introduces an uncertainty in the background predic- 
tion that varies between 24% and 34%, and an uncer- 
tainty on the mSUGRA signal efficiencies between 15% 
and 17%. Uncertainties related to the modeling of the 
initial- and final-state soft gluon radiation in the sim- 
ulated samples translate into a 3% to 6% uncertainty 
on the mSUGRA signal efficiency, and uncertainties on 
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the background predictions that vary between 8% and 
10%. An additional 10% uncertainty on the diboson and 
top quark contributions accounts for the uncertainty on 
the predicted cross sections at NLO. A 2% uncertainty 
on the measured Drell-Yan cross sections, relevant for 
Z/7*+jets and IF+jets processes, is also included. The 
total systematic uncertainty on the SM predictions varies 
between 31% and 35% as the jet multiplicity increases. 
Various sources of uncertainty in the mSUGRA cross sec- 
tions at NLO, as determined using prospino, are consid- 
ered. The uncertainty on the PDFs varies between 10% 
and 20%, depending on the mSUGRA point considered. 
Variations of the renormalization and factorization scales 
by a factor of two change the theoretical cross sections 
by 20% to 25%. 

Figure [T] shows the measured Ht and $ T distributions 
compared to the SM predictions after all final selection 
criteria are applied. For illustrative purposes, the Fig- 
ure indicates the impact of a given mSUGRA scenario. 
The measured distributions are in good agreement with 
the SM predictions in each of the three final states con- 
sidered. In Table II, the observed number of events and 
the SM predictions are presented for each final state. A 
global x 2 test, including correlations between systematic 
uncertainties, gives a 94% probability. 



Events in data (2 fb 



> 2 jets > 3 jets 
18 38 



> 4 jets 
45 



SM predictions 



4.4 ±2.0 13.3 ±4.6 
1.3 ±1.2 7.6 ±4.1 



QCD jets 
top 

Z^uP+jets 3.9 ±0.9 5.4 ±1.4 
Z/7* ->Z+Z-+jets 0.1 ±0.1 0.2 ±0.1 
W^lv+]ets 6.1 ±2.2 10.7 ±3.1 
WW,ZW,ZZ 0.2 ±0.2 0.3 ±0.2 



15.3 ±7.1 
22.1 ± 7.0 
2.7 ±0.7 
0.1 ±0.1 
7.7 ±2.2 
0.5 ±0.2 



total SM 



16 ± 5 37 ± 12 48 ± 17 



TABLE II: Number of events in data for each final state com- 
pared to SM predictions, including statistical and systematic 
uncertainties summed in quadrature. 

The results are translated into 95% C.L. upper limits 
on the cross section for squark and gluino production in 
different regions of the squark-gluino mass plane, using 
a Bayesian approach [29| and including statistical and 
systematic uncertainties. For the latter, correlations be- 
tween systematic uncertainties on signal efficiencies and 
background predictions are taken into account, and an 
additional 6% uncertainty on the total luminosity is in- 
cluded. For each mSUGRA point considered, observed 
and expected limits are computed separately for each of 
the three analyses, and the one with the best expected 
limit is adopted as the nominal result. Figure [2] shows the 
observed and expected 95% C.L. upper limits as a func- 
tion of squark and gluino masses, compared to mSUGRA 
predictions, in four regions within the squark-gluino mass 



DATA (L = 2.0 fb' 1 ) mSUGRA signal: 
SM M s = 350 GeV/c 2 

Total Syst. Uncertainty M = 385 GeV/c 2 
SM + mSUGRA A.= 0, tan(5 = 5, n<0 




;10' 



&10 



H x > 280 GeV 
Er> 90 GeV 
> 4 jets 




H T > 330 GeV 
120 GeV 
> 3 jets 




H T > 330 GeV 
> 180 GeV 
> 2 jets 
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FIG. 1: Measured Ht and $ T distributions (black dots) 
in events with at least two (bottom), three (middle), and 
four (top) jets in the final state compared to the SM predic- 
tions (solid lines) and the SM±mSUGRA predictions (dashed 
lines). The shaded bands show the total systematic uncer- 
tainty on the SM predictions. 



plane. Cross sections in the range between 0.1 pb and 
1 pb are excluded by this analysis, depending on the 
masses considered. The observed numbers of events in 
data are also translated into 95% C.L. upper limits for 
squark and gluino masses, for which the uncertainties on 
the theoretical cross sections are included in the limit cal- 
culation, and where the three analyses are combined in a 
similar way as for the cross section limits. Figure [3] shows 
the excluded region in the squark-gluino mass plane. For 
the mSUGRA scenario considered, all squark masses are 
excluded for M~ < 280 GeV/c 2 , while for M~ = M~ 

r2 



<i a 

masses up to 392 GeV/c 2 are excluded. Finally, for 
M~ < 400 GeV/c 2 gluinos with M~ < 340 GeV/c 2 are 
excluded. This analysis extends the previous Run I limits 
from the Tevatron by 80 GeV/c 2 to 140 GeV/c 2 . 

In summary, we report results on an inclusive search 
for squarks and gluinos in pp collisions at ^fs — 1.96 TeV 
in events with large $ T and multiple jets in the final 
states, based on 2 fb -1 of CDF Run II data. The mea- 
surements are in good agreement with SM predictions for 
backgrounds. The results are translated into 95% C.L. 
upper limits on production cross sections and squark and 
gluino masses in a given mSUGRA scenario, which sig- 
nificantly extend Run I results. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 



7 



Academy of Finland. 



w 2 

,-, 10 

°3 1 

gio 1 



o 10 
t/3 



o 

Cio- 



ur 



-Theoretical uncen a i nlics no 
in ihc calculation of ihc lim 


included 




- - NLO Prospino 

syst. uncert. (PDF © Ren.) 

expecled limif 95% C.L. 

— observed limit 95% C.L. 


" M g = 230 GeV/c 2 






= 370 GeVAr 








A = 0, tanp = 5, \i<0 








<^ L = 2.0fb" 1 


_ M g = M q 






= 460GeV/c 2 



300 400 500 200 300 400 500 
- 2l M„ [GeV/c 2 ] 



M q [GeV/c 2 ] 



FIG. 2: Observed (solid lines) and expected (dashed lines) 
95% C.L. upper limits on the inclusive squark and gluino 
production cross sections as a function of M~ (left) and 
M~ (right) in different regions of the squark-gluino mass 

plane, compared to NLO mSUGRA predictions (dashed- 
dotted lines). The shaded bands denote the total uncertainty 
on the theory. 
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FIG. 3: Exclusion plane at 95 % C.L. as a function of squark 
and gluino masses in an mSUGRA scenario with Ao = 0, 
fi < and tan/3 = 5. The observed (solid line) and expected 
(dashed line) upper limits are compared to previous results 
from SPS H0| and LEP [H| experiments at CERN (shaded 
bands) , and from the Run I at the Tevatron [j| (dashed-dotted 
line) . The hatched area indicates the region in the plane with 
no mSUGRA solution. 
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